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We report the results of certain integrations of quantum-theoretic interest, relying, in this regard, 
upon recently developed parameterizations of Boya et al of the n x n density matrices, in terms of 
squared components of the unit (n — l)-sphere and the n x n unitary matrices. Firstly, we express 
the normalized volume elements of the Bures (minimal monotone) metric for n = 2 and 3, obtaining 
thereby "Bures prior probability distributions" over the two- and three-state systems. Then, as a 
first step in extending these results to n > 3, we determine that the "Hall normalization constant" 
(Cn) for the marginal Bures prior probablity distribution over the (n — l)-dimensional simplex of 
the n eigenvalues of the n x n density matrices is, for n = A, equal to 71680/7r^. Since we also 
find that C3 = 35/n, it follows that C4 is simply equal to 2^^C3/tv. (C2 itself is known to equal 
2/n.) The constant C5 is also found. It too is associated with a remarkably simple decompositon, 
involving the product of the eight consecutive prime numbers from 3 to 23. We also preliminarily 
investigate several cases n > 5, with the use of quasi-Monte Carlo integration. We hope that the 
various analyses reported will prove useful in deriving a general formula (which evidence suggests 
will involve the Bernoulli numbers) for the Hall normalization constant for arbitrary n. This would 
have diverse applications, including quantum inference and universal quantum coding. 

PACS Numbers 03.65. Bz, 03.67.Hk, 02.70.Lq, 02.60.Jh 
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I. INTRODUCTION 



We make use of recently proposed parameterizations (cf. |^) — in terms of squared components of the unit 
(n — l)-sphere and the n x n unitary matrices — of the n x n density matrices. First (in sec. ||), we derive (prior) 
probabiUty distributions of particular interest over both the three-dimensional convex set of two-state quantum systems 
and the eight-dimensional convex set of three-state quantum systems. These distributions are the normalized volume 
elements of the corresponding Bures metrics on these systems. Hall |^ (cf. has contended that such distributions 

correspond to "minimal-knowledge" ensembles, that is the most random ensembles of possible states. In particular, for 
the two-dimensional quantum systems, he argues that the Bures metric provides such a minimal-knowledge ensemble, 
since it "corresponds to the surface of a unit four-ball, i. e., to the maximally symmetric space of positive curvature 
. . . This space is homogeneous and isotropic, and hence the Bures metric does not distinguish a preferred location or 
direction in the space of density operators" Somewhat contrarily though. Slater 0| has reported results (based 
on the concept of comparative noninformativities of priors, first expounded in ||^) that indicate the Bures metric 
generates ensembles that are less noninformative than other (monotone) metrics of interest. 

The Bures metric fulfills the role of the minimal monotone metric ||9|"|ll|] , and has been the focus of a considerable 
number of studies |^,0,|l2|-^ . "An infinitesimal statistical distance has to be monotone under stochastic mappings" ||, 
p. 786]. All stochastically monotone Riemannian metrics are characterized by means of operator monotone functions. 
Among all (suitably normalized) operator montone functions f{t) with /(I) = 1 and f{t) = tfit~^\ there is a minimal 
and a maximal one p^ . (The concept of a minimal metric was apparently introduced by Zolotarev in his extensive 
paper, "Metric distances in spaces of random variables and of their distributions" sec. 1.4], but it is not entirely 
clear that the meaning there is the same as in the terminology "minimal monotone metric".) 

We should bear in mind, though — as emphasized by Petz and Sudar ]^ — that, in strong contrast to the classical 
situation, in the quantum domain there is not a unique monotone metric, but rather a (nondenumerable) multiplicity 
of them. Their comparative properties need to be evaluated, before deciding which specific one to employ for a 
particular application. We repeat the concluding remarks of Petz and Sudar: "Therefore, more than one privileged 
metric shows up in quantum mechanics. The exact clarification of this point requires and is worth further studies." 

We have previously reported (in terms of parameterizations other than that of Boya et al Q) the Bures 
probability distribution for the two-state systems, and also for an imbedding of these systems into a four- 
dimensional convex set of three-state systems, but the result below ([III ) for the full eight-dimensional convex set of 
three-state systems is clearly novel in nature. In fact, in ]^, sec. II. E] we discussed certain (unsuccessful) efforts in 
these directions (although the volume element of the maximal monotone metric — which is not strictly normalizable 
— proved more amenable to analysis there). 

In sec. in, we determine certain necessary elements for extending the work reported in sec. |^ to the higher- 
dimensional quantum systems (n > 3). This involves finding the normalization constant (C„), explicitly first discussed 
by Hall eq. (25)], for the marginal Bures prior probability distribution over the (n — l)-dimensional simplex of the 
n eigenvalues of the n x n density matrices. These constants are found to exhibit quite remarkable number-theoretic 
properties. It would, therefore, certainly be of substantial interest to find a general formula for C„. Knowledge of the 
value of Cn, together with that of the invariant Haar element for SU (n) — apparently presently available, however, in 
suitably parameterized form (cf. j2^j2^) for n < 3 25| — would allow one to construct the Bures prior probability 
distribution itself for the n-level quantum systems. 

In an extensive study [ p6[ , Krattenthaler and Slater examined (in the framework of the two-staie systems) the 
hypothesis that the normalized volume element of the Bures metric would function in the quantum domain in a 
role parallel to that fulfilled classically by the "Jeffreys' prior" — that is, the normalized volume element of the 
unique monotone/Fisher information metric [ pTp^ ]. In particular, they were interested in |26| in the possibility of 
extending certain (classical) seminal results of Clarke and Barron |2^j3^. They did conclude, however, contrary to 
their working hypothesis, that the normalized volume element of the Bures metric does not in fact strictly fuUfill the 
same role as the Jeffreys' prior (in yielding both the asymptotic minimax and maximin redundancies for universal 
coding/data compression), but it appears to come remarkably close to doing so (cf. Fig. In sec. [V, for the cases 
n = 2 and 3, the "quasi-Bures" prior probability distributions are presented that appear to fulfill this distinguished 
information-theoretic role. 



II. BURES PROBABILITY DISTRIBUTIONS OVER THE N x N DENSITY MATRICES 

Boya et al ^ have recently "shown that the mixed state density matrices for n-state systems can be parameterized 
in terms of squared components of an (n — l)-sphere and unitary matrices" . The mixed state density matrix (p) is 



2 



represented in the form, 



p = UDU\ (1) 

where U denotes an SU{n) matrix, its conjugate transpose and D a diagonal density operator, the diagonal entries 
(di's) of which — being the eigenvalues of p — are the squared components of the (n — l)-sphere. Thus, for n ~ 2, 

--{-^T Jen) (^) 

and, for n = 3, 

'cos2 0/2 sin^ 61/2 \ 

D=\ sin^ 0/2 sin^ 61/2 < 61, < tt. (3) 

cos^ 6'/2 / 



(Note the differences in the ranges of angles used in the two cases. This will be commented upon in sec. [II B ) 
Biedenharn and Louck have presented ||3l|, eq. (2.40)] the parameterization of an element of SU{2), 

in terms of the Pauli matrices (di's) and three Euler angles — < a < 27r, </3<7r,0<7<27r — with an 
associated invariant Haar measure eq. (3.134)], 

dl72 = idad7sin/3d/3. (5) 
8 

Byrd (cf. |^^) has extended this approach to SU{3). He obtains 

U{a, (3, 7, K, a, b, c, C) = e'^'°'e'^^'^e'^'^e'^''^e'^''''e'^^^e'^"'e'^<"^ , (6) 

where Ai denotes one of the eight 3x3 Gell-Mann matrices |3^]. The corresponding invariant element is 

dfis — sin 2/3 sin 26 sin 2Ksin^ KdQ!d/?d7dKdad6dcdC, (7) 

with the eight Euler angles having the ranges, 

< a,7,a,c < TT, < ^9, 5, k < 7r/2, < ( < \/3. (8) 

For our purposes, the Euler angle 7 for the case n = 2 and c and ( in the case 71 = 3 are irrelevant, as they "drop 
out" in the formation of the product (0). (I thank M. Byrd for this important observation.) So, we will employ below 
the appropriate conditional versions of these invariant measures (^ and (0) — the condition (a technical statistical 
term) corresponding, of course, to the ignoring of the indicated angles. 
The Bures metric itself is expressible in the form eq. (10)] 

where \i > denotes the eigenvectors of the n x n density matrix p, < j|, the corresponding complex conjugate (dual) 
vectors, and the d's are the associated eigenvalues. The parameterization of Boya et a/is, then, particularly convenient, 
since the eigenvalues and eigenvectors of p are immediately available. Our chief concern must, then, be to compute the 
complete Jacobian of the transformation to the set of parameters of Boya et al. (We shall note for further reference 
the occurrence of the term di + dj in (^. This, is of course, simply proportional to the arithmetic mean, (di + dj)/2. 
By replacing this term by (twice) the exponential/identric mean (p3| ) of di and dj, that is 2I(d,i, dj), we shall olatain 



the particular "quasi-Bures" distributions described in sec. IV 
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A. The Bures case n = 2 



For the case n — 2, the volume element of the Bures metric is proportional to the product of the inverse of the 
square root of the determinant of p (or, equivalently, the determinant of D) with two Jacobians. The first Jacobian (in 
line with the familiar practice in the theory of random matrices eq. (3.3.5)]) is itself the product of (di — ^2)^ and 
the (conditional) invariant element (^). The second Jacobian, sin0/2, corresponds simply to the transformation from 
cartesian coordinates to the squared polar coordinates employed in (|^). Simplifying and normalizing the full product, 
we arrive at the probability density for the normalized volume element of the Bures metric over the three-dimensional 
convex set of two-state quantum systems. This density is 

cos'^ sin 3 tt 
dpBures-2{0, «, P) = 5 — -d0dad/3 O<0<-, 0<a<27r, 0</3<7r. (10) 

TT^ 2 

The expected values of the eigenvalues are, then, 1/2 ± A/Stt. 



B. The Bures case n = 3 

For the case n ~ 3, the volume element of the Bures metric is equal to the product of: (i) two Jacobians again, 
one of which now has the form ((di — d2){di — d^){d2 — ds)) multiplied by the (conditional) invariant measure (7), 
while the other, (cos | sin'^ | sin(/))/2, corresponds to the transformation to squared spherical coordinates used in (S); 
and (ii) the reciprocal of the product of the square root of the determinant of p (or, equivalently, of D) and the 
difference between the sum of the three principal minors of order two of p (or, equivalently, of D) and the determinant 
itself. Since = (cos | sin^ | sin0)/2, it can be seen that considerable cancellation occurs between the numerator 

and the denominator of the full product. The normalization of the resultant volume element required considerable 
manipulations using MATHEMATICA (basically involving reducing the problem to the simplest possible form at each 
stage of the integration process). We obtained the following Bures prior probability density over the eight-dimensional 
convex set of three-state (spin-1) quantum systems, 

35u 

1287r4(35 + 28 cos 6* + cos 261 - 8 cos 2(/)sin'' 

where 



dpBuresA^^, a, (3, J, K, a, b) = , , — n 4 es, d6ld0dr^3, (11) 



u — sin - ((35 + 60 cos e + 33 cos 29) cos ^ - 8 cos Scpsin^ -f), (12) 

and the conditional invariant element (cf. (0)) is 

dfl^ = sin 2/3 sin 2b sin 2Ksin^ KdQ!d/3d7dKdad5. (13) 

The eight variables have the previously indicated ranges ((^), In Fig. |l|, we display the two-dimensional marginal 
probability distribution of (|l^) over the parameters 9 and cj) (which are invariant under unitary transformations of p). 




. 6 

prob. densit^Q _ 4 
0.2 

f 

^1 

2' 

3^0 

FIG. 1. Bivariate marginal Bures prior probability distribution over the variables d and <j> for three-state systems 
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Let us note that the fully mixed state — corresponding to the 3x3 diagonal density matrix with entries equal to 
1/3 — is obtained at (j) = 7r/2, = 2 cos~-^(l/\/3) ~ 1.91063. The probabihty density ( pl| ) is zero at this distinguished 
point, as well as along the loci 9 = and (fi = tt/2. (Wherever at least two of the eigenvalues of p or, equivalently D, 
are equal, the density is zero.) 

The one-dimensional marginal probability density (Fig. |^), obtained by integrating (^ij) over all variables except 6, 

is 

35 6 30 f^O 6 

dpBuresa(0) = — ("1533 + 2816 COS - - 1988 cos 6i + 1152 cos — - 447cos20 + 128 cos — ) sin^ -AO. (14) 
256 2 2 2 2 



prob. density 
. 1 



. 08 



.06 



. 04 



. 02 




0.5 1 1.5 2 2.5 3 
FIG. 2. Univariate marginal Bures prior probability distribution over the variable 9 for three-state systems 



The relative maxima of this density are located at .914793, 2.2795 and tt, while the relative minima are at 0, 1.59995 
and 2.61732. 

The one-dimensional marginal probability density (Fig. ^ , obtained by integrating ( |ll] ) over all variables except 0, 

is 

d^B^resisW = 3^(cot0cscV(llOlOO48Otan"i [cot ^Jcos^^ ^ - 26880(792(27r - 0) cos (/) -I- 87r (15) 

(55 cos 30-^3 cos 5(/)) -f (/)(495 cos 20- 220 cos 30 + 66 cos 40 - 12 cos 50 + cos 60) ) + 16885656 sin 20 
-^5069937 sin 40 -f 167012 sin 60 - 3(41395200 -I- 124 sin 80-4 sin 100 + sin 120)) ) d0. 



prob. density 
1 . 5 
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-0 . 5 
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0.5 1 1.5 2 2.5 3 



FIG. 3. Univariate marginal Bures prior probability distribution over the variable for three-state systems 
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The limits of ( p^ as 4> approaches and tt are both equal to |^ « .707355. (The density is highly oscillatory in 
the vicinity of these boundary points.) MATHEMATICA does, in fact, perform a symbolic/exact integration of jl^), 
yielding the result 1. (Note that for our analyses below for n > 3, we have found it necessary to rely upon numerical 
integrations, although the results obtained do appear to indicate that exact solutions exist, which in principle might 
be found with a powerful enough computer.) However, several warning messages are generated en route to this 
result, concerning indeterminate expressions and inconsistencies in the arguments of MeijerG functions (which are 
very general forms of hypergeometric functions) 



III. HALL NORMALIZATION CONSTANTS FOR MARGINAL BURES PROBABILITY DISTRIBUTIONS 
OVER THE (TV l)-DIMENSIONAL SIMPLEX OF THE TV EIGENVALUES OF THE TV x TV DENSITY 

MATRICES 

Let us note that Hall eq. (24)] has, in fact, given an explicit formula for the volume element of the Bures metric 
on the n X n density matrices. This is (converting to the notation used above), 

where the real part of the ij-entiy of the diagonalizing unitary matrix U in formula (^ is represented by Xij and the 
imaginary part by iiyij. Since in the parameterization of Boya et al [Q, which we have employed, one uses not these 
x's and y's, but rather the Euler angles parameterizing the unitary matrix U, we have been compelled to replace the 
differential elements dxijdyij in (^|) by the corresponding conditional form (dfi) of the invariant (Haar) measures (|^) 
and (0). One can, then, confirm that our presentation and results are fully consistent with the use of (16), bearing in 



mind the unit trace requirement that di + . . . + d„ = 1. Hall eq. (25)] also expressed the marginal Bures probability 
distribution over the space of n eigenvalues of p as 

PBures-.nidl, ■ • • , d„) - ^» (di...d„)V2 U d, + d, ' 

We shall report the values of the "Hall constants" C„ for n = 2, 3, 4 and 5, immediately below. (The results for n > 2 
are, apparently, new.) 



A. The Hall constants for n = 2 and n = 3 



If, for consistency with our further results for n > 2, we take 9 € [0, tt] and not [0, n/2] as in Q], then we find that 
C2 equals 2/7r w .63662 (cf. |3[ eq. (30)]). From the results of the analysis in sec. p B| , we are able to determine, for 
the first time, apparently, that C3 — 35/tt « 11.1408. Let us note here that 35 is, of course, simply the product of 
proximate or neighboring prime numbers, that is, 35 = 5 • 7. 



B. The Hall constant for n = 4 and associated methodology 



To continue the full line of research reported here for the cases n = 2, 3, to n > 3, it would be useful to extend 
the work of Byrd and Sudarshan p^ , p5t o n the Euler angle parameterization of SU{3) to such higher n. However, 
computation of the Hall constants C„ (|17|) does not depend on parameterizations of SU{n). We have, in fact, been 
able to obtain exceedingly strong numerical evidence that C4 is, in fact, equal to 71680/7r^ « 7262.7. Let us, first, 
mention some methodological considerations useful in deriving this result (and, in general, C„, n > 4). 

The parameterizations of Boya et al Q of the diagonal 2x2 and 3x3 matrices differ, in that in the 2x2 case (|^) 
only matrices in which the (l,l)-entry is at least as great as the (2,2)-entry are generated (due to the restriction of the 
angular parameter, 9, to the range [0, 7r/2]), while in the 3x3 case no order is imposed on the diagonal entries (the 
parameters 9 and 4> both varying freely between and tt). Now, in performing (the apparently necessary) numerical 
(as opposed to symbolic) integrations to obtain the Hall constants C„ (n > 3), it seems to be considerably more 
computationally effective to integrate over only those diagonal matrices in which (say) the (l,l)-entry is no less than 
the (2,2)-entry, which in turn is no less than the (3,3)-entry, etc. (This helps to minimize troublesome oscillations.) 
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Then, the result can be multiphed by the number (n!) of permutations of n objects to yield C„, since the result of the 
integration must be invariant under any other of the n\ — 1 possible orderings (permutations) that can be imposed 
on the diagonal entries of the n x n diagonal matrices. In precisely this manner, we were able to obtain (using the 
numerical integration of interpolating function command of MATHEMATICA) the result 71680. OOOOOI/tt^ for C4. 
This we take as overwhelming evidence that C4, in fact, equals 71680/7r^, particularly so, since 71,680 has the highly 
structured prime decomposition of 2^^ • 5 • 7. (It is also interesting to note that 35 = 5 • 7 appears in the numerator 
of C3, that is, 35/n, or, to the same effect, C4 = 2^^C3/7r. 

To illustrate the procedure followed, let us first parameterize the 4x4 nonnegative diagonal matrices of trace unity 
in the following fashion (cf. (||), (||)), 

\ 

I sin2 6'/2cos2 0/2 ° 1 nx^ 

sin^ 61/2 sin^ 0/2 cos2C/2 ' 

V sin2 6'/2sin2 0/2sin2C/2/ 

Then, the (truncated) region of integration employed above — corresponding to successively nonincreasing diagonal 
entries — can be described as C G [0,7r/2],(/) e [O,/(^)],0 e [Q,f{(j))], where f{x) = 2cot~^ (cosa;/2), rather than 
C'pj^^ G [Oi'''']i in the apparently suggested parameterization of Boya et al [y, which would yield all possible 
diagonal matrices, without regard to the ordering of their elements. 

If we construct a similar truncated region of integration in the case n = 3, then we find that the expected values of 
the eigenvalues are .802393, .181878 and .0157299. In the n = 2 case, the analogous values are (as previously noted), 
l/2±4/37r, that is .924413, and .0755868. 



C. The Hall constant for n = 5 



We have also attempted to compute C5, in the manner of sec. [II B, with the use of MATHEMATICA. We obtained 
(using the Gauss-Kronrod integration method with a working precision of twenty-one digits, rather the machine 
precision of sixteen) the result 2342475135. OO/tt^. Now, it is most interesting to note (particularly, in light of our 
results for C„, n < 5, that 

2, 342, 475, 135 = 2in^^2^, = 21(3 • 5 • 7 • 11 • 13 • 17 • 19 • 23), (19) 

where Pi denotes the i-th prime number (taking the sequence of primes to be 2,3,5, . . .). Thus, we have acquired 
strong evidence that, in fact, C5 = 2342475135/7r^. (As a simple exercise, we looked at the one hundred thousand 
consecutive integers containing 2,342,475,135 as their midpoint, and computed all their prime decompositions. All 
the others had at least one prime factor greater than 23.) The seemingly independent factor of 21 in (|9|) will also 
apparently be found below in the (odd) case n = 7. 



1. Prime factorials 

In the product of the primes less than or equal to p is denoted p#. (The issue there, as in several of the 
works cited there, was to test p# ± 1 for primality (cf. @).) Let us point out a 1952 article entitled "Tables of 
logarithms of the prime factorials from 2 to 10007" (a synopsis of which can be found in Mathematical Reviews 16, 
112f). It is noted that by the prime number theorem the ratio of the sum of the logarithms of the primes from 2 
to p to p itself approaches 1 as p ^ 00. So, if it eventuates that the general formula for C„, at least for odd n, contains 
a term of the form p(n)#, where p{n) is a prime as a function of n, which grows indefinitely large with n itself, then 
it should be possible to asymptotically replace p(n)# by e^*^"^. ("A version of the prime number theorem states that 
the product of the primes less than x is asymptotically [citing the well-known treatise Theorem 434]], but the 
error term is notoriously large, so it is probably unrealistic to expect to be able to compute far enough to get within 
the necessary epsilon" Q.) We also observe that with the use of Wilson's theorem Q, {p — 1)! = — Imodp, one 
could express p^ in terms of the (more) standard factorial function p\ . (In , the primes are defined in terms of 
factorials.) 

Let us point out that EUinas and Floratos have recently studied the "prime decomposition" of an n x n density 



matrix into a sum of separable density matrices with dimensions determined by the coprime factors of n 43 1 (cf . 
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D. Preliminary investigations of the Hall constant for n > 5, with the use of quasi-Monte Carlo integration 



1. n = 6 

Of course, as the dimensionality of the n-state quantum systems increases (that is, n itself increases), the numerical 
integrations required to sufhciently narrow estimates of the corresponding Hall constant become increasingly more 
difficult, and it is hard to judge what is precisely the optimum numerical/programming strategy to employ. Following 



the methodology outlined in sec. Ill B , based on the ordering of the eigenvalues, MATHEMATICA did yield (using the 
standard default options) an estimate of Cq, representable in the form, 1.4616286 • IO^^/tt'^ (although diagnostics as to 
inadequate precision were issued during the course of the computation). This result, coupled with our observation of 
the pattern of C„ for n < 6, might lead us to speculate that the numerator of Cg is either the scventccn-digit number, 
14, 616, 907, 579, 654, 144 = 2'^'^ ■ 17^ • 23 or 14, 623, 504, 649, 420, 800 = 2'*2 . 5^ • 7 • 19, with the denominators, in both 
cases, being n^. 

When we employed the quasi-Monte Carlo (Halton-Hammersley-Wozniakowski) procedure chap. 3] of 
MATHEMATICA, to numerically integrate over a hypercube of volume tt^ (corresponding, thus, now to no particular 
distinguished ordering of the six eigenvalues), we obtained (with no diagnostics at all being generated in two separate 
analyses — having set maxima of ten and fifty million sample points) a result of the form 1.536355674- IO^^/tt"^. To 
a very high accuracy, this numerator can be approximated by 15, 363, 556, 773, 986, 304 — 2^^ • 3** • 13^ • 37 • 113 for the 
numerator of the presumptive value of Cg (with denominator, again, tt^). Sloan and Wozniakowski have noted that 
recently "Quasi-Monte Carlo algorithms have been successfully used for multivariate integration of high dimension d, 
and were significantly more efficient than Monte Carlo algorithms" However, in comparing the two sets of results 
here for n = 6 it is very important to bear in mind that MATHEMATICA sets its precision and accuracy objectives 
much lower when Monte Carlo procedures are employed j4^]. (Of course, these default values can be reset, but in the 
preliminary analyses reported here, they have not been, though we intend to do so in future studies.) 

2. n = 7 

The quasi-Monte Carlo MATHEMATICA procedure produced an estimate (with no accompanying diagnostics, 
having set a maximum of four million sample points) of C7 as 2.4811899 • The numerator of this fraction 

(cf. (|l9|)) could be approximated (to a relative error of less than three-tenths of one percent) by the product of 21 
and the (eighteen) consecutive primes from 2 to 61 (in the notation of this is 61#), that is, C7 — 21Ill^^Pi/TT^ . 
However, this degree of accuracy does not at all seem satisfactory, particularly in light of the proximate results for 
n — 6 and 8, though odd n's appear to present greater computational challenges. 

5. n = 8 

The quasi-Monte Carlo procedure estimated Cg as 4. 1836777028 -IO^^/tt^ (though, unlike the n = 6, 7 cases, a failure 
to converge was reported — that is, with the preassigned use of at most ten million sample points). Nonetheless, 
this outcome can be fit to a very high accuracy (less than one part in one hundred million) by taking Cs to be 
289.32-5-132-31-47-61/7r''. 

E. n > 9 

For the case n = 9, convergence was not obtained with the use of six million sample points. The result given 
was 7.631832917- IO'^^/tt'*. We have also made tentative attempts, using the quasi-Monte Carlo procedure again, 
to estimate C„ for n = 10 and 11 (but it appears that considerable investment of computer resources is needed for 
sufficiently satisfactory answers). For n — 10, based on a maximum of two million sample points, convergence was 
not obtained and the result 6.334733996 • IO^^/tt^ reported. (No related decomposition was immediately apparent.) 
For n = 11 and 12, using a maximum of one million sample points in both cases, the results 3.254198489 • IO^^/tt^ 
and 5.218327334- IO^'^/tt^ were gotten (without convergence or any obvious associated simple prime decompositions, 
however) . Our impression is that the computations are considerably more difficult for the odd values of n than for the 
even, which may be some reflection of the simpler formulas displayed above for even n (cf. Pq,E9[). (For n > 7, the 
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MATHEMATICA compiler was unable to handle the small numbers appearing in the calculation, and then proceeded 
with the use of the uncompiled evaluation, leading to slower running times p^.) Of course, there exists a wide 
range of possible approaches to numerical integration problems of this kind, including, certainly, the use of alternative 
programming languages, in particular, FORTRAN. The trade-offs between these various options need to be assessed. 

As a reference point, against which one can attempt to compare the (reciprocals of the) several values of C„ above, 
let us recall that the area of an rt-sphere of radius r is given by 27r"/^r"~^/r(rt/2) |Q. Also of similar interest is 
Euler's formula 

C(27i) = i— f— (2^)2"B2„, 71-1,2,3,... (20) 



2(2n; 



where C(s) is the Riemann zeta function, X)n>i (necessarily rational) n-th Bernoulli number 

pl| . Vol. I, pp. 75 and 211]. (The values of ({n) for positive odd integers n, however, have not been expressed in 
such a simple form |5^, Vol. Ill, p. 1695]. Infinite series representations are, in fact, reported in p8| , p9[ .) Initial 
attempts to find an explanatory formula for the sequence of (integral) numerators of C„, using on-line programs of 
C. Krattenthaler ("Rate") App. A] and of N. Sloane ( "superseeker" ) Q, did not succeed. However, we did 
eventually find a somewhat intriguing connection (further buttressed by some related analyses, discussed below in 
sec. IIIG) between the results here and sequence A035077 of ||5^, which gives the denominators of partial sums of 
B2n- The numerator of: (1) C2 is (rather trivially) twice the first entry in this sequence; (2) C3 is one-half the fourth 
entry; (3) d is 2^° = 1024 times the fourth entry; and (4) C5 is 63 times the thirteenth entry of AG35G77. (We 
note that Donaldson has found simple proofs of various formulas for symplectic volumes involving Bernoulli numbers 



F. Average von Neumann entropy of 71-state quantum systems with respect to Bures prior probability 

distributions 

As one application of these computations of the Hall constants (C„), let us note that with respect to the Bures 
probability distribution ( p^ the average von Neumann entropy, — Trp log p, is exactly 2 log 2 — 7/6 .219628 nats for 
the two-state systems and, now using numerical integration, .507937 nats for the three-state systems (cf. p8| ). (Since 
we employ the natural logarithm here, the unit of information is the nat, which is equivalent to l/log2 w 1.4427 
bits.) This latter result, to a high degree of precision — that is, to ten significant places, can, in fact, be written 
as 3 log 3 — 3917/1405. We also computed the Bures average entropy for the four-state quantum systems, obtaining 
.751771 nats. To eight significant places, this can be written as 4 log 4 — 32135/6704. Similarly, for n — 5, using our 
knowledge of C5, we obtain an average entropy of .954103 nats. This is closely approximated (to at least nine places, 
according to our calculations) by 5 log 5 — 40045/5648. 



G. Auxiliary analyses of variations of Hall integrals and the role of Bernouilli numbers 

Since it appeared to be quite challenging to determine the Hall constants (C„) for n > 5, we thought that it might 
be revealing (possibly helpful in deriving a general formula for arbitrary n), as well as being of independent interest, 
to investigate more tractable variations. To do so, we replaced the exponent two i n (|l6| ) and (|^ by either one 
(corresponding to real quantum systems) or four (for quaternionic quantum systems) |}3|] , and other positive integer 
values of less immediate physical interest, as well. 

For the n — 3 (spin-1) quaternionic case (that is, using an exponent of four), the counterpart of the previously 
derived (sec. Ill A ) Hall constant 35/7r is 1616615/2267r = 5 • 7 • 11 • 13 • 17 • 19/2267r. (The numerator here is one-half 
the tenth member of the sequence A035077, comprised of denominators of partial sums of the Bernoulli numbers i?2n 
Hi.) When we use an exponent of six, the result is 100280245065/882527r = 3 • 5 • 7 • 11 • 13 • 17 • 19 • 23 • 29 • 31/882527r. 
(This numerator is precisely the seventeenth member of A035077.) 

If we use an exponent of one (corresponding to the case of real quantum systems), the result is l/47r, and if we 
employ an exponent of three, the normalization constant is 105/1287r = 3 • 5 • 7/1287r. (The numerator here is the fifth 
member of A035077.) For an exponent of five, we have 15015/81927r = 3 • 5 • 7 • 11 • 13/81927r, the numerator being 
three times the seventh entry of Sloane's sequence. For an exponent of seven, the result was the eleventh member of 
the sequence A035077, that is, 969969, divided by 2621447r. 

When an inquiry was made of Neil Sloane as to whether to his knowledge there were any published discussions of 
this sequence, he replied "No, I was just looking at various sequences of important rationals, and thought that the pair 
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A035078/A035077 should be in the database". (However, he later pointed out that the von Staudt-Clausen Theorem 
p9| , p. 10] was relevant to questions involving sums of Bernoulli numbers.) For the case n = 3, MATHEMATICA 
rejected our efforts to compute any further exact integrals having integer exponents greater than seven. (It would 
appear that the use of an exponent of eight would be associated with the octonionic quantum systems [pO[.) 

For the analogous set of variations with n = 2, use of odd exponents in (|l6|) and ( p7| ) lead to divergent results. For 
an exponent of two, the result is 27r, for four, 8/37r, for six, 16/57r, for eight, 128/357r, .... 

For comparable scenarios based on n = 4, we were unable to proceed with exact integrations. Our numerical 
computation of the analog of the Hall constant employing an exponent of unity (the real quantum case), yielded 
TT^/. 079271. But we were unable to determine if this result bore any relation to the sequence A035077. 



IV. QUASI-BVK'ES PROBABILITY DISTRIBUTIONS OVER THE N xN DENSITY MATRICES 

In line with the work reported in [p6| , it would be of interest to obtain formulas for the averages over the eight- 
dimensional convex set of 3 x 3 density matrices (p) with respect to the Bures prior probability distribution (^Tj) of 
the m-fold tensor products of p. As m ^ oo, the relative entropy of these products with respect to the averaged 
3™ X 3"* density matrix gives us the (Bures) asymptotic redundancy for the universal quantum coding of three-state 
systems. 



A. The quasi-Bures case n — 2 

In |26[ | and further yet unreported work, the Bures prior probability density (|l^) for the two-state systems was 
found to closely resemble the (what we term "quasi-Bures" ) probability density, 

dpgnasi:2(6',a,^) = •226231(tan'''='^''6'/2)cos6'cot6'sin/3d6'dad/3, O<0<^, 0<a<27r, < (3 < tt (21) 

which yields, it appears, both the asymptotic minimax and maximin redundancies (as "Jeffreys' prior" [ pTf does 
classically [ p9|j30| ]). This common value, if one ignores the error term, as appears to be legitimate, is (31ogm)/2 — 
1.77062, while the Bures probability distribution ( p^ ) has been shown to be, incorporating the error term, associated 
with an asymptotic redundancy of (31ogm)/2 — 1.77421 -I- 0{l/^/rri) [^, p. 29]. In general, for any probability 
distribution w{6),6 e [0,7r/2], the asymptotic redundancy for the two-state quantum systems takes the form, 

3 TTl 1 

- log 2 log sin 6 + 2(scc 9) log tan log w{6) + o(l). (22) 

2 27r 2 2 

Standard variational arguments can, then, be used to show (ignoring the error term, the legitimacy of which seems 
plausible, but has not yet been rigorously justified) that the particular w{9) yielding both the maximin and minimax 
redundancies is simply proportional to the ^-dependent part ot pquasi:2iS , a, (3) , indicated in (^l|). 

The reciprocal of the corresponding "Morozova-Chentsov function" c{x, y) PJlC|] for ( |2l] ) is the exponential or 
idcntric mean [ |6l] , eq. (1.3)] of x and y, 

I{x,v)^e-\x^/yy)^y, x^y (23) 

{I{x,x) — x), while for the Bures (minimal monotone) metric, it is the (more commonly encountered) arithmetic 
mean [x -f ?;)/2. The associated operator monotone functions P,|l0|] are f{t) = (1 -I- t)/2 for the Bures metric, 
and /(t) = the metric giving (^l|). (The Morozova-Chentsov functions fulfill the relation c{x,y) = 

l/yf{x/y).) Perhaps the exponential mean arises in this context because the von Neumann entropy is the logarithmic 
relative entropy, and of course the exponential and logarithmic functions are inverses of one another. This leads 
us to speculate that if one were to employ, following |l^], the "quadratic relative entropy" or the "Bures relative 
entropy" instead, then, in the parallel universal coding context, the minimax/maximin would be achieved by the 
means corresponding to the new forms of inverse functions. While the logarithmic relative entropy is based on the 
operator convex function, g{t) — — logt, the quadratic form relies upon g{t) = (t — 1)^ and the Bures form on 
g{t) ^{t- l)V(t + 1) 0. Here, f{t) = (i - l)V(g(i) + tg{l/t)). 

In Fig. ^, we jointly display the univariate marginal probability distributions of (^0|) and (^), revealing that they 
closely resemble one another, with the quasi-Bures distribution assigning relatively greater probability to the states 
more pure in character {9 < .443978). 
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prob. density 
1 . 4 




0.25 0.5 0.75 1 1.25 1.5 
FIG. 4. One-dimensional marginals of the [n — 2) Bures probability distribution and the quasi-Bures distribution jil] ) 
(assuming greater values for 9 G [0, .443978]), which yields both the asymptotic minimax and maximin redundancies for the 
universal quantum coding of two-state systems 



B. The quasi-Bures case n = 3 

The three-state counterpart of (|2l|) (that is, the probability distribution associated with the exponential/identric 
mean, rather than the arithmetic mean, as for such Bures distributions) is 



dpquas^■.3i^, 4>, 7, &) = .000063495u(tan(i+"°= '^V/2)csc'' -csc^ 



(24) 



(cos — tan-) "+'" (sm — tan-) d^dtpdUs, 



where u is given in (|T^), df^a in ( [T^ ) and 

i; = 2 + 6cos6', w = cos (61 - (/)) - 2 cos (/) -I- cos (6* + 0). (25) 
The corresponding two-dimensional marginal probability distribution over the variables 9 and (f> is exhibited in Fig. |^. 



prob. densit^-'^ 
0.2 




3 

FIG. 5. Bivariate marginal of the quasi-Bures probability distribution (^^ over the variables and for three-state systems 

As would be anticipated from Fig. this figure closely resembles Fig. |l|. In Fig. ^, we show the result obtained by 
subtracting the bivariate marginal Bures probability distribution shown in Fig. |l] from its quasi-Bures counterpart in 
Fig. |. 



11 




FIG. 6. DifFerence obtained by subtracting the three-state Bures bivariate marginal probabihty density in Fig. |l] from the 
quasi-Bures one, associated with the exponential/identric mean, displayed in Fig. H 



C. Quasi-Bures counterparts of Hall (Bures) normalization constants 



If we replace the term dj + dk, occurring in the denominator of the expression (|l7|), by (twice) the ex- 
ponential/identric mean ( p3| ) of dj and dk, that is, I{di,dj), the resultant expression becomes a formula for 
Pquasi-.nidi, . . . , dn) , now interpreting C„ to be the normalization constant for the corresponding quasi-Bures proba- 
bility distribution. For n = 2 then (taking 9 G [0,7r] and not [0,7r/2] as in ||l|), rather than 2/7r w .63662, we obtain 
.769427, and for n = 3, instead of 35/7r .089759, we find .138681. 



V. CONCLUDING REMARKS 

We hypothesize that the (full eight-dimensional) quasi-Bures probability distribution ( p4[ ) associated with Fig. ^ will 
furnish the common asymptotic (minimax and maximin) redundancies for universal quantum coding in that higher- 
dimensional setting, paralleling the result (not yet fully formally demonstrated, however) for the three-dimensional 
convex set of two-state systems (cf. ||62|). In this regard, it might prove computationally convenient, as a heuristic 
device, to replace the quasi-Bures probablity distributions by their (closely approximating) Bures analogues, since 
certain exact (symbolic) integrations are achievable (at least, for the cases n = 2, 3) with the Bures distributions, but 
apparently not with the quasi-Bures ones, for which numerical methods seem to be necessary. 

Although a parameterization of SU(4) is not relevant, as already noted, to the computat ion o f the Hall constant C4, 



that is, 71680/7r^, and to that of the corresponding average (von Neumann) entropy (sec. Ill F ), it would be essential 
in investigating the universal coding of four-state quantum systems, since the tensor products of 4 x 4 density matrices 
have to be calculated and it would, therefore, be necessary to implement formula (^. In a personal communication, M. 
Byrd has indicated that he has undertaken the (challenging) task of developing such an (Euler angle) parameterization 
of SUiA). 

Such a parameterization of SU (4) — in conjunction with the knowledge, acquired here, of C4 — might also prove 
of value in estimating the volume of 2 x 2 separable quantum states and lead to numerically more stable results 

than those reported in Table 1] , since the "over-parameterizations" of the unitary matrices used there could then 
be avoided, due to the "dropping out" (as pointed out for n = 2, 3 immediately after (H) above) of certain Euler 
angles in the formation of the product (^. We also note that in |6j] (cf. (6^]), the Bures (minimal monotone) metric 
was found to yield higher a priori probabilities of entanglement than other monotone metrics (in particular, the 
Kubo-Mori-Bogoliubov and maximal ones). Presumably, even in any computationally improved form of analysis, this 
conclusion would be unaltered. 

We have also pursued a traditional (pseudo-random number) Monte Carlo approach to estimating the Hall normal- 
ization constants for 2 < n < 16. However, the degrees of precision attained were not satisfactory. (For discussions of 



the comparative computational complexities of the pseudo- and quasi-Monte Carlo methods, see 47|66|].) 

We are presently attempting to obtain a more precise estimate of the Hall constant Cg, in particular, using non- 
Monte Carlo (that is, adaptive) integration methods. Our best current estimate of Cq is 1.534836628- IO^^/tt'^. The 
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numerator can be very well approximated by 15348366279966720 = 2^"^ • 3 • 5 • 7^ • 11 • 13 • 17, which seems more 
satisfactory than the results reported in sec. IIIDl, based on the quasi-Monte Carlo procedure. Now, it is most 
interesting to note that this numerator is precisely 2^^^ • 7 times the ninth entry of the sequence A035077, we have 
repeatedly referenced above. We are also compelled to observe that our educated conjecture as to the numerator of 



(sec. [II D 2) is exactly forty- two times the thirty-third entry 
published list does not extend this far — of A035077. 



which we had to compute ourselves, since Sloane's 
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